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Cooling of neutrons stars proceeds mainly via neutrino emission; As an example we 
study the modified neutrino bremsstrahlung process nn — » nnvv. The radiated en- 
ergy is small compared to other scales in the system. Hence one can use low-energy 
theorems to compute the neutrino emissivity in terms of the non-radiative process, 
ie the on-shell T matrix. We find that the use of the T— matrix as compared to 
previous estimates based upon one pion exchange leads a substantial reduction of 
the predicted emissivity. 



1 Introduction 

The thermal evolution of neutron stars is dominated by the weak interaction 
and in particular neutrino interactions with the hadronic matter. The hadronic 
information is contained in the socalled current-current correlator II (q), where 
q = (q,u>) is the four-momentum of the neutrinos. In general one can distin- 
guish two different regimes: 

• neutrino scattering (space-like, ui < q) 

• neutrino-pair emission (time-like, lu > q). 

In the latter case the quasi-particle response vanishes, i.e., the one-body pro- 
cess n ^ n + v + D is forbidden by energy-momentum conservation, and the 
socalled URCA process 

• n—>p + e~+D 

p + e~ — ► n + 

is forbidden unless the n,p and e~ fermi momenta satisfy the unequality 
p p F + p e F > p F , which is very unlikely. 

Therefore the bremsstrahlung process can only take place in the presence of 
spectator nucleon (referred to as modified processes), e.g.: 

• n + N — > n + N + v + v (N = n,p) 

• n + n — > n + p + e~ + D (+ inverse) 

In the pioneering work of Friman and Maxwell □ these processes were com- 
puted in the extreme soft neutrino limit in Born approximation with the NN 
interaction presentated by a Landau type interaction plus a one-pion exchange 
part. 

With respect to this approach several questions can be asked; how accurate is 
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Figure 1: Leading order diagrams for soft bremsstrahlung 



the use of Born approximation, what is the contribution of other mesons, and 
how important are relativistic effects? 

ft is the aim of the present study to address these issues by computing the 
emissivity using a low-energy theorem on the basis of the observation that 
the neutrino radiation is very soft compared to other scales in the system 
(u> fts T = 1 MeV); this allows one to use an empirical T— matrix, fully deter- 
mined by phase shifts. 

2 Soft electro-weak bremsstrahlung 

First we consider the "Soft-photon" amplitude in free space 
The original soft-photon theorem states that the first two terms in an expansion 
of the electromagnetic bremsstrahlungs amplitude in photon four-momentum 
q are fully determined by the amplitudes of the non-radiative process, u 

M = A/lj + B + 0(w) 

This result can be generalized in several ways, e.g. to virtual bremsstrahlung, 
(q 2 > 0) and also to the case of the weak axial vector current. Here we shall 
restrict ourselves to the leading order, the A term, in which case only radiation 
off external legs contributes (see Fig. 1), and the amplitude can de expressed 
as 

M„ = TS{p! - q)Y li + r„S(pi + q)T + {1^2} 

The weak vector and axial vector vertices Y are given by low-energy neutral 
current Hamiltonian, 

with the hadron current given by = Yli=n,p V'*7 ,i (^w,i — TbCa,*)^*' Hence 
in the non-relativistic limit 

(vector) « g v -f^ -> S^q 
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T^(axial) fa 5^757^ -> 3^ct (i = p, n; g\ = -g\) 

The two Feynman propagators (corresponding to prior and post emission) are 
given by 

1 ±m ±1 1 

S(p ±q) = —— « = — - « ±-(1 + 0(p.q/mu,)). 

r ){p±q) — m p.q ui — p.qjm uj 

Hence in this order the amplitude M can be expressed as a commutator of T 
with the space component of the axial weak current operator 

M i = ^[T,S i ] 

UJ 

where S — (<ji + o"2)/2, and in this order there is no vector contribution. 
Which terms survive the commutator [T, r wcak ] ? Looking at the structure 
of the NN amplitude given in the appendix one sees that there will be non- 
vanishing contributions from tensor, spin-orbit and quadratic spin-orbit terms 
in T. Thus one may conclude that only nn spin correlations contribute to the 
emissivity. However, we note that in the past the effective nn interactions were 
restricted to purely local interactions and hence the spin-orbit interactions were 
not considered. As an illustation we compute the cross section for NN electro- 
weak bremsstrahlung in free space as a function of the relative nn momentum 
(in the nn cm system). In fig. 2a we show the results for full T— matrix and 
its separate components. In fig. 2b the cross section for the full T— matrix is 
compared with the result for 7r, tt + p and also a meson exchange. 
We note that the OPE exchange (which forms the basis of most standard 
"cooling scenarios" overestimates the results obtained with the full T— matrix 
by about a factor 5 for relative momenta in the range appropriate for neutron 
matter at normal density (j>f = 1.3fm _1 ) (a similar result has been found in 
□). It is also seen that the contribution from the OPE tensor force is largely 
canceled by that from p exchange, and that the sum of OPE and ORE is much 
closer to the T -matrix result than OPE . 



3 The np case 

The np case is slightly different, since the axial vector couplings of neutron and 
proton have opposite values g p A = —g\ = —1.25. Hence in this case the net 
result can be expressed as the sum of a commutator and an anticommutator, 
where the latter corresponds to the exchange diagrams 

Mi = i^L([T dir : Di] + {T exch , Di}) 
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Figure 2: Cross section for weak nn bremsstrahlung for lu= lMeV as a function of the relative 
momentum p = (pi — p2)/2; top: full T-matrix (solid line), tensor+spin-spin (dashed), spin- 
orbit (dotted) and quadratic spin-orbit (dashed-dottcd line); bottom: OPE (solid), OPRE 
(dotted), OPR.SE (long dashed), direct OPE (short dashed), direct OPRE (dashed-doubly 
dotted) and full T (dashed-dottcd line). 
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where D = (<ti — (t 2 )/2, and T dlr ' exch are denned in the appendix. 
In this case also the socalled Landau terms in the T-matrix, of the form 
<?CTi-C2and g' 'cr1.a2Ti.T2 contribute to the anti-commutator. In practice the 
proton fraction in a neutron star is quite small and therefore we find that the 
relative contribution from np is quite unimportant. 

4 Emissivity in medium 

The emissivity in the medium with given density and temperature T can be 
computed from the amplitude in two ways: (i) via the (imaginary part) of 
the current-current correlator, or (ii) using the Fermi golden rule. The former 
approach is more general, but the latter is simpler to use and is sufficient in the 
present case. In lowest order in the virial expansion it amounts to convoluting 
the free space matrix elements \M 2 \ with the hadronic Fermi-Dirac functions 

/eft eft 
rf 3 Pi — — (wi + lu 2 )8\ Pi +p 2 -P3-P4-qi~ q 2 )F\ < Mj > | 2 , 

where F — /i/2(l— /a)(l— fi)- In practice the integrals are simplified by taking 
the absolute values of the nucleon momenta equal to the fermi momenta. If 
we consider the ratio of the emissivity calculated with the full T matrix and 
the one from the OPE£ we obtain about a factor 4-5 reduction. This is inline 
with the conclusion inta for pure neutron matter. 

Other medium that need to be considered in further work are 1) replacement 
of the free T- by a G-matrix, which takes into account the Pauli blocking and 
other medium effects, 2) medium renormalization of the axial coupling vertex, 
3) inclusion of higher order medium effects such as dressing Green functions. In 
the latter case one has to take care to conserve the symmetries (CVC, PCAC) 
of the problem (which are conserved in the present case) . 
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Appendix: Structure of the on-shell NN amplitude 

The on-shell T-matrix is fully determined by the nn phaseshifts. First we 
consider nn; in a covariant approach one has 

T on = T a (s,t, u)u(p 2 )tt a u(p 2 ) x u(pi)n a u(pi) 
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with 

tt a = {l,ev,„, 7^,75,7^75} 

Non-relativistically 

T = T C + T s ai.a 2 + T T Si 2 + T S oL.S + T Q Q 12 

In case of np there are twice as many components which can be distinguished 
by isospin / = 0, 1 , i.e. T dir = (T° + T 1 )/2 and T exch = (T° - T 1 )/2. 
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